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<N ■ ABSTRACT 

Q^ ■ We have computed theoretical models of circumstellar disks for the classical 

Be stars k Dra, j3 Psc, and v Cyg. Models were co n struct ed using a non-LTE 



% 



^ ■ radiative transfer code developed by ISigut fc JonesI (120071 ) which incorporates 

I . a number of improvements over previous treatments of the disk thermal struc- 

5_i ■ ture, including a realistic chemical composition. Our models are constrained by 

^ . direct comparison with long baseline optical interferometric observations of the 

Ha emitting regions and by contemporaneous Ha line profiles. Detailed com- 
parisons of our predictions with Ha interferometry and spectroscopy place very 
j^ . tight constraints on the density distributions for these circumstellar disks. 
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1. Introduction 

Classical Be (or B-emission) stars are rapidly rotating, hot stars with optical spectra 
that show hydrogen emission lines, and, fr equently , emiss ion lines from singly ionized metals. 



It has been recognized since the days of IStruvd (jl93ll ) that recombination in a flattened 



disk of circumstellar gas can reproduce the basic features of the spectroscopically observed 
Ha line profiles of Be stars. However, it is now generally accepted that this model is too 
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simplisti c. The detailed mechanisni(s) that cr eates and maintains Be star disks remains 
unclear ( Porter fc Rivinud 12003 : Owockil 2003), although rapid rotation of the central B 
star certainly plays a role (JTownsend et al.l |2004| ) . In addition to this lack of a successful 
dynamical model, the evolutionary status of Be stars and related stars (such as B[e] stars) 
is not well understood, and possible evolutionary connections between these groups are 
ambig uous. These problems are compounded by the addi tional uncertainty in their rotation 



rates (jCranmerll2005l : iTownsend et al.ll2004j : IOwockj|2003f) and lack o f und erstanding of how 



rotation and evolution interplay. The suggestion by ITownsend et al.l (120041 ) that the rotation 
rates of Be stars may be systematically underestimated, combined with new interferometric 
data, has rekindled interest in finding suitable dynamical models for these stars. 

Until recently, the bulk of the observational evidence for the presence of disks surround- 
ing Be stars has been spectroscopic and polarimetric in nature. The spectroscopic evidence 
includes the profiles of the emission lines observed in the optical and IR spectra, as well as 
continuum excesses at infrared wavelengths. The net linear polarization observed in Be stars 
arises via electron scatterin g in the non-spherical distribution of gas, as first suggested by 
Coyne fc Kruszewskil (1l969l ). However, interferometric observations that spa tially resolve the 



disks in several wavelength regimes are becoming available, (see for example I Chesneau et al. 



20051 : iMcAlister et al.l l2005l : iTycner et al.ll2006l ). Interferometric observations, in concert 



with other observables, allow the disk physical properties to be determined with greater 
accuracy and this may contribute to the development of more successful dynamical models. 



Tycner et al.l (120061 ) collected long-baseline, optical, interferometric observations of the 
classical Be stars 7 Cas and Per at the Navy Prototype Optical Interferometer (NPOI) 
and interpreted these observations using simple models of the intensity distribution on the 
sky, such as uniform disks, r ings, or disks with a Gaussian intensity distribution. In this 
paper, we extend the work of iTycner et al.l (120061 ) by computing detailed theoretical models 
for the intensity distribution on the sky produced by the central B star and disk of a Be 
system. We then compare these models to observations collected at the NPOI. 



The theoretical models are computed with the bedisk code developed by lSigut &: Jones 



(I2OO7I ). Assuming a density distribution of the disk gas, the code can compute the tempera- 
tures in the disk given the energy input from the central stars' photoionizing radiation field, 
by enforcing radiative equilibrium. As the disk density model contains several adjustable 
parameters, values for these parameters can then be extracted from the match to obser- 
vations. We present new observations and detailed models for the Be stars k Dra, (3 Psc, 
and V Cyg. We also use near-contemporaneous observations of the Ha spectral lines for 
these stars as additional constraints on the density model. As we shall demonstrate, the 
additional spectroscopic constraint of detailed line profiles is often an important ingredient 



-3- 



in selecting among models consistent with the interferometric observations. The overall goal 
of this study is to place tight constraints on physical conditions in the circumstellar regions 
of these Be stars. 



2. Theory 

The theoretical disk n aodels presented in th is investigation were constructed using the 
BEDISK code developed by lSigut &: Jonej (120071 ). A brief overview of the cod e and assump- 
tions relevant to this work are presented below. The reader is referred to ISigut fc Jones 
( 120071 ) for more details. 



The circumstellar disk is assumed to be axi-symmetric about the star's rotation axis and 
symmetric about the mid-plane of the disk. In this case, cylindrical geometry is appropriate 
and we use R for the radial distance from the star's rotation axis and Z for the height above 
the equatorial plane. The radial density distribution in the equato r ial plane is assu r ned t o 
be given by an - R~" p ower-law following the works of I Waters! (Il986l ): ICote fc Waters! ( 11987! ): 
Waters fc Cote! (11987! ). At each radial distance from the rotation axis, the vertical structure 
of the disk is determined by the requirement of hydrostatic equilibrium in the Z direction 
which balances the Z component of the star's gravitational acceleration with the gradient of 
the gas pressure. 

Typical values for the radial power-law index n for theoretical models of Be stars usuall y 
fall in the range of 2 to 3.5 based on models that fit the IR continuum (jWatersI !l986! ). 
Models with the highest values of n have a faster decrease in the density distribution with 
increasing R. The density distribution within the disk is determined by an assumed value 
of the density at the stellar surface in the equatorial plane, po, by the power-law index, n, 
and by the assumption of pressure support perpendicular to the equatorial plane. Since po 
and n significantly affect the theoretical predictions, we conduct a 2-dimensional parameter 
search by varying these input parameters over all reasonable values for comparison with 
interferometric and spectroscopic observations. 

Given this density model, the bedisk code is used to find the temperature structure 
of the disk by enforcing radiative equilibrium. The microscopic heating and cooling rates 
implied by a gas with a solar chemical composition are computed and the temperature that 
balances heating and cooling is found. To compute the atomic level populations required 
by this procedure, the statistical equilibrium equations are solved for each chemical element 
included, accounting for the bound-bound and bound-free coUisional and radiative processes 
that set the rates in and out of each atomic level. 
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To compute the hydrogen hne profiles, we solved the transfer equation along a series of 
rays through the star plus disk system as viewed at an angle i relative to the line-of-sight 
(where z = 0° is pole-on, and i = 90° is equator-on). The disk is assumed to be in pure 
Keplerian rotation, and the equatorial rotational velocity of the star was chosen so that the 
measured v sin i of the star was recovered for the adopted value of i. For rays terminating on 
the stellar surface, we adopted a photosphe ric Ha line profile computed using the synthe 
code and the Stark Broadening routines of iBarklem fc Piskunovl (120031 ): these profiles are 
based on the hydrogen pop ulations c ompu ted for the appropriate LTE, line-blanketed model 
atmosphere adopted from iKuruczl (jl993l ). For rays passing through the disk, the equa- 
ti on of radiative tr a nsfer was solved along the ray using the short-characteristics method 
of lOlson &: Kunaszl (jl987l ). To compute the Ha opacity and emissivity, the hydrogen level 
populations computed for the thermal solution were used. Here, hydrogen was represented 
by a 15 level atom and all implied radiative and coUisional bound-bound and bound-free 
processes from and between the 15 levels were included. The radiative bound-bound rates 
were trea ted by the esc ape probability approximation in which the net radiative bracket for 
the line (lMihalaslll978l ) is replaced by a single-flight escape probability. While the escape 
probability approximation can provide a reasonable description of the line thermalization, 
the use of escape probabilities is an important approximation. The bedisk code assumes 
static escape probabilities based on the compl ete redistribution oy er a Doppler profile for 
all lines (for a justification of this procedure see lSigut &: JonesI 120071 ). If the functional form 
of the escape probability function is changed, for example, to a simple r~^ dependence on 
optical depth for r >■ 1, the predicted equivalent width of Ha can change by ~ 2 A (due to 
the change in the thermal structure of the disk and in the hydrogen level populations). Thus 
in matching to the observed Ha profiles, the largest uncertainty is likely from the uncertainty 
in the theoretical profile and not from the observational uncertainties (see Section 3 for an 
estimate of the observational uncertainties). 

Our approach allows the full dependence of the Ha line emissivity and opacity on the 
physical conditions t h rough out the disk to be retained. The Stark broadening routines of 
Barklem &: Piskunovl (120031 ) were also used to compute the local Ha line profile throughout 
the disk. This is important as these routines can handle high disk densities for which the 
assumption of a simple Gaussian or Voigt profile for the Ha line would be invalid. The final 
Ha line profile from the unresolved system was obtained by summing over the rays weighted 
by their projected area on the sky. The Ha line profile was computed over a wavelength 
bin of ±20A from line centre at each area element on the sky. To obtain the final Ha 
line profile to compare with spectroscopic observations, the profile was convolved with a 
Gaussian of FWHM of 0.656 A to bring the resolving power of the computed profile down to 
10^ to match the observations. The predicted interferometric visibilities also followed from 
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the same numerical model. The monochromatic Ha image of the system projected on the 
plane of the sky was computed by integrating the stellar continuum and Ha line flux over a 
150 A spectral window centered at Ha to represent the interferometric spectral channels (see 
Section 3). Thus the calculation of the Ha profiles and the Ha interferometric visibilities 
used the same model calculations. 



3. Observations 

We obtained long-baseline interferometric and spectroscopic observations of three Be 
stars, K Dra, (3 Psc, and v Cyg. The interferometric observations were acquired using 
the Navy Prototype Optical Interferometer (NPOI) utili zing baselines ranging f rom 18.9 



to 64.4 m in length. The NPOI is described in detail in lArmstrong et al.l (119981 ) and the 



technique used to extract interferometric obs ervables from t he sp ectral channel containing 



the Ha emission line has been discussed by iTycner et al.l (120031 ). Table [T] shows all the 



interferometric observations obtained for the three stars and the number of data points 
obtained on each night in the spectral channel containing the Ha line. The interferometric 
observations of all stars were obtained over a period of less than a month, with (3 Psc and 
V Cyg observed during 2005, and k Dra observed during 2006. 

To complement our interferometric observations we have also obtained near contempo- 
raneous spectroscopic observations of all three targets. The spectroscopic observations were 
obtained using an Echelle spectrograph at the Lowell Observatory's John. S. Hall telescope 
and covered the spectral region around the Ha line. The properties of the spectrograph 
and the descriptio n of the reduction methods have been described elsewhere (see § 3.2 in 



Tycner et al.ll2006l . and references therein). The continuum signal-to- noise ratio (SNR) near 
Ha is more than 200 in the spectra of all three stars. However, the equivalent width and 
the shape of the line is determined with respect to the normalized continuum level. Be- 
cause the line profile must be normalized by a smoothly varying function that is fit to the 
continuum (or where one thinks the continuum is located), the choice of this normalizing 
function affects the shape of the Ha line. In addition, the presence of telluric lines, as well 
as their temporal variability from night to night (and season to season) can affect how well 
the continuum level is determined; we expect this to be the dominant source of uncertainty 
associated with the shape of the observed Ha line. We have estimated the magnitude of this 
effect by analyzing a set of reduced spectra all obtained from the same raw spectrum, but 
all obtained with slightly different quadratic functions fitted to slightly different continuum 
regions. The variations in the resulting normalized continuum were at the 2 to 3% level, and 
to be conservative we adopt a 3% uncertainty for our continuum level determination and, in 
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turn, the shape and equivalent width measure of the observed Ha emission hne. 

For (3 Psc and v Cyg the spectra were obtained on 2005 Sep 16, and for n Dra the 
spectrum was obtained on 2006 Mar 17, and therefore in all three cases the spectroscopic 
observations overlapped the interferometric observing runs (Tabled]). The equivalent widths 
of the observed Ha spectra lines are -15.3 A, -24.8 A, and -22.2 A for f3 Psc, v Cyg, and 
K Dra, respectively. (We adopt the usual convention that a negative in the quoted value of 
the equivalent width indicates a net line flux above the continuum.) 



4. Results 

There is considerable uncertainty in the assigned spectral types for Be stars due to 
rapid rotation and spectral variability. The effects of gravity darkening, and the potential 
obscuration of portions of the stellar surface by the disk, further compound the uncertainty 
in classification and the corresponding stellar parameters. Since the stellar parameters fix 
the photoionizing radiation field that is assumed to be the sole source of energy input into the 
circumstellar disk, realistic stellar parameters are essential to the construction of reasonable 
disk models. For each of the stars in this study, we searched the literature in order to find the 
most reasonable stellar parameters available. See Table [2] for the adopted stellar parameters 
for K, Dra, j3 Psc, and v Cyg. The rationale for adopting the particular stellar parameters is 
explained and compared with other published values in each subsection for a given star. 



4.1. K Draconis 



K Dra is a known spectroscopic bi nary and gian t star of intermediate spectral type. It is 
variab le on a varie ty of timescales; see iHiratal (119951 ) for a detailed discussion of these varia- 



tions. iSaadI (120041 ) modeled the central star of n Dra by comparing a variety of spectroscopic 
and photometric observations to a grid of non-LTE models. We have adopted their derived 
best-fit stellar parameters for our study. We us ed their parameters of Te//, logg, and mass 
to calculate the stellar radius. iGies et al.l (120071 ) also used these parameters for the interpre- 
tation of -ft^^-b and observat ions of k Dra collected at the CHARA Array interferometer. The 
parameters of ISaadI (120041 ) are also in reas onable agreement with stellar parameters listed 
for the spectral class of B5III in ICoxl (J2000l ). 



Figure [T] shows the density parameter grid for the disk of k, Dra. This grid was selected 
after running a series of models over a much wider range in po and n on a coarser grid. 
A comparison of the observed and predicted Ha equivalent widths allowed us to determine 
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Tabl e 1: Observing Log for Interferometric Observations 
UT Date 

2005 Aug 28 

2005 Aug 29 

2005 Aug 30 

2005 Aug 31 

2005 Sep 1 

2005 Sep 5 

2005 Sep 13 

2005 Sep 14 

2005 Sep 15 

2005 Sep 16 

2005 Sep 17 

2005 Sep 26 

2006 Feb 25 

2006 Feb 26 

2006 Mar 5 

2006 Mar 6 

2006 Mar 14 

2006 Mar 16 

2006 Mar 19 

Total: 



K Dra f3 Psc v Cyg 


^ of points) (# of 


Doints) (t^ of points 




12 


S 




16 12 




20 20 




20 20 




20 20 




( 


B 


4 




16 16 




23 20 




20 20 




21 24 




( 


B 


3 




16 28 


30 






6 






52 






24 






36 






76 
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Table 2: Adopted stellar parameters 



Star 


HR 


log^ 


Spectral 


Stellar Radius 


Stellar Mass 


^e// 


vsini 


Reference 


Name 


Number 




Type 


R0 


Mo 


K 


km/s 




K Dra 


HR4787 


3.5 


B6IIIpe 


6.4 


4.8 


14000 


170 


1 


/3Psc 


HR8773 


4.0 


B6Ve 


3.6 


4.7 


15500 


90±15 


2 


V Cyg 


HR8146 


4.0 


B2Vne 


4.7 


6.8 


19800 


173±10 


3 



References. — fl lSaadI (|2004l) ; (;2 lLevenhageii fc Leisteii (|2004l) : (3) iNeined (|2005l ) 
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Fig. 1. — The density parameter grid for k Dra shown as a function of po (g/cni^) and n. 
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Fig. 2. — The x^ values for computed models as a function of n and po (g/cm^) for k Dra. The 
contour plot shows the reduced x^ values in steps of 1.0. Only a region of reduced x^ < 6 
is shown. The darkest regions on the figure represent models with reduced x^ < 2. The 
four white plus signs correspond to models with the lowest values of reduced x^ determined 
from interferometry. The diamonds correspond to models that have predicted Ha equivalent 
widths ± 2 A of the observation. 
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Fig. 3. — The predicted Ha equivalent widths as a function of log po, and the power-law 
index, n, for the entire grid of models. The lightest coloured symbols are closest to the 
value of the observed equivalent width of -22.2 A . The models that predict Ha equivalent 
width within ± 2 A of the observation are represented by white circles. The models that 
correspond to predicted Ha emission greater than the observation are indicated by the plus 
sign within the symbol. 
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Fig. 4. — The observed Ha line for k Dra compared with 4 models with n and po as indi- 
cated in the panels. These models correspond the best-fit from both interferometry and Ha 
modeling. The model profiles are the solid lines and the circles represent the observed line. 
All model profiles were obtained with a disk inclination of 35°. The reduced x^ values from 
interferometry for these models are 1.18, 1.29, 1.18, and 1.18, from left to right respectively. 
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Fig. 5. — The interferometric data and model with n = 4.2, and po = 1.5 x 10"^° g/cm^ for 
K, Dra. The dashed hne represents the central star and the solid line is obtained by taking a 
Fourier transform of a synthetic image of the preferred model. 
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a finer grid to sample for detailed comparisons with interferometric observations. The fi- 
nal density parameter grid consisted of 278 models. We have compared the output from 
each model to the interferometric observations using the technique described in detail by 



Tycner et al.l (120081 ). The reduced x^ values corresponding to these models are presented in 
Figure [2] as a function of po (g/cm'^) and n. The darkest patches in this Figure represent 
the models with reduced y^ < 2. The white plus signs correspond to models with the lowest 
values (ranging from 1.18 to 1.29) of the interferometric reduced x^. In total, there are 25 
models with reduced x^ <2 for this star, therefore, it is necessary to further constrain our 
models in order to determine the best-fit model. 

Figure [3] shows the predicted Ha equivalent widths as a function of log po, and the 
power-law index, ra, for the entire grid of models. The predicted Ha profiles were constructed 
assuming an inclination of i = 35°. We note that this value of i is consistent with other 
values deter mined from sp ectroscopy and polarization analyzes presented in the literature. 



For examp l e. IClarkd (119901 ). estimates i = 23° based on the observed polarization angle, and 
Juza et al.l (jl99ll ) find i = 35° — 45° based on spectroscopy. We also experimented by varying 
i and found that an inclination angle for this system of 35° seemed to produce profiles that 
matched the observed Ha line. We also note that although changing the value of i by ± 10° 
changes the shape of the Ha profile, the effect on the equivalent width is < 1 A. The lightest 
coloured symbols are closest to the value of the observed equivalent width of -22.2 A. There 
are 14 models that predict an Ha equivalent width within ± 2 A of the observation and 
these are represented by white circles. Models which correspond to predicted Ha equivalent 
widths greater than the observed emission have a plus sign within the corresponding symbol 
in Figure [31 

We note that this subset of 14 models does not represent the 14 best models using 
interferometry alone. Although this set matches the observed Ha equivalent width within ± 
2 A, using interferometric observations as a constraint this same set corresponds to reduced 
X^ from 1.18 (the minimum value) up to 16.0. In fact, 6 of the 14 models in this set had a 
reduced x^ greater than 2.0. It is clear from a comparison of Figure [2] and Figure [3] that the 
results from Ha interferometry and Ha line profile are complimentary but that each provides 
extra information to help reduce the degeneracy in finding the best-fit model. 

Multiple models, of pairs of n and po, that produce Ha profiles with equivalent widths 
that match the observations within uncertainties can be found. However, not all of these 
are the best-fits in terms of interferometry. That is, it is possible to match the equivalent 
width for a range of n by an appropriate adjustment in po, but that these combinations do 
not necessary have the correct density distribution as a function of radial distance from the 
star. However, by combining results from spectroscopy and interferometry, we can find the 
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best pairs of n and po for k Dra. 

Next, we compared the subset of 8 models which represent the best-fits from both 
interferometry and Ha equivalent width. A subset (4 models with the lowest reduced x^ 
from interferometry) of the corresponding Ha profiles is shown in Figure HI Notice the 
variation in the shape of the profile as n is increased from 2.7 in the left-most panel to 4.2 
in the right-most panel. The predicted equivalent widths are -21.9, -23.5, -22.4, and -20.3 
A, with the corresponding x^ values from interferometry of 1.18, 1.29, 1.18, and 1.18. The 
model with n = 2.1 does not have the the correct line profile shape compared to the observed 
profile. There is too much absorption at line centre, and the profile is clearly too narrow in 
the wings. All of the predicted profiles in Figure H] have an absorption feature at line centre. 
This feature is the most pronounced for low values of n and gradually decreases for models 
with higher values of n. If we increase the inclination, this absorption feature becomes more 
and more pronounced as the system is viewed closer and closer to edge on. If we decrease 
i, so the system is viewed close to pole-on, the absorption feature remains for the models 
with the lowest values of n. In addition, with n = 2.7, the wings of the predicted profile are 
not broad enough for any reasonable values of i. The panel with the model corresponding 
to n = 3.8 and n = 4.0 in Figure H] shows that the predicted line is too strong and this 
results in an Ha equivalent width which is too large. Models corresponding to n > 4.2 have 
profiles with wings wider than the observation. The model that corresponds to the profile 
with n = 4.2 matches the shape of the observed profile best and has a minimum reduced 
X^ of 1.18 from interferometry. We prefer this model with n = 4.2 and p„ = 1.5 x 10^^° 
g/cm^ in terms of both spectroscopy and interferometry. Figure O shows the interferometric 
observations compared with this model for k, Dra. 

We note that as previously discussed the uncertainty due to the approximations in the 
computational code (see Section 2) are larger than those in the observed Ha profiles (see 
Section 3). We searched the literature to find the best possible stellar parameters for the 
disk systems studied in this investigation (see Section 4.1, 4.2, 4.3). However, generally 
for Be stars, not only are the stellar parameters not well established but for some stars, 
even the spectral type may not be accurate. We now wish to test the dependence of model 
parameters on our results. We experimented by varying the stellar parameters in turn, to 
see how these changes affect the predicted Ha equivalent widths and the interferometric 
fits for the star k Dra. Similar types of uncertainties would be expected for the other two 
systems presented in this paper. One might expect that the Tg// of the central star, which 
supplies the disk with energy, may be the source of the greatest uncertainty. We changed 
the Tgff by ± 1000 K and reproduced all of the profiles. This resulted in a change in the 
equivalent width of the Ha line by a maximum of ~ 2 A. To assess the effect of uncertainties 
in the fundamental stellar parameters combined, we performed a Monte Carlo simulation 
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in which new stellar parameters were randomly realized. These parameters were used to 
compute a new circumstellar disk model, corresponding to our preferred model of n = 4.2 
and po = 1.5 x 10~^° g/cm^. This procedure was repeated to generate new Ha profiles and 
equivalent widths. Errors assumed for the stellar parameters were ± 25% in mass, ± 25% 
in radius, and ±10.7% in Te// (corresponding to ±1500 K); each stellar parameter was 
assumed equally probable within its errors. Using a Gaussian distribution for the errors, we 
found our results did not change significantly. The mean of Ha equivalent width was —20 A 
with a standard deviation of 2.65 A. 



We have one star in common with the work of iGies et al.l (120071 ): k Dra. They use a disk 
density model which is very similar to ours and derive two sets of parameters (n, po) based 
on fits to K'-band interferometric visibilities: n = 0.16 ± 0.43, po = 2.8 x 10 
single star model, and n = 0.67±0.36, po = 6.2 x 10 



13 gcm_Mor a 



13 g ^T^ 3 |-Qj, ^ binary model. IGies et al. 



(120071 ) prefer the binary fit to the single star model because of the peculiar solution they 
obtain in the single star case; namely a very small exponent in the radial density power law 
(n = 0.16) and a best solution with i = 72° ± 18° which is at odds with the resolved Ha line 
profile. Nevertheless, the particular properties of the introduced binary are not constrained 
by any independent observations. Both of these solutions seem at odds with our results. 
It is well docum ented that the Ha equival e nt wid t h and the infrared excess are correlated 
(see for example Ivan Kerkwijk et al.l Il995l ) . ISaadI (120041 ) present a figure whi ch sh o ws th e 
change in Ha equivalent width from the early 1970's to 2004 for n Dra (see ISaadI (120041 ) . 
figure 4(d)). Although there is certainly scatter of data values in their figure, one can easily 
see from this figure that the Ha equivalent width varies by ~ 15 A (from about - 5 to 
about -20 A) in a period of ~ 6 years. This is a substantial change in the Ha emission 
line and we expect there would be a correspondi ng increase in th e infrared excess. The 
interferometric observations for k Dra presented in iGies et al.l (120071 ) were obtained in April 
(3 dates) and December (2 dates) of 2005. Their spectroscopic observations were obtained 
in April 2005. Our conteniporane ous observations were obtained in 2006. The difference in 
our results and iGies et al.l (120071 ) is partially accounted for by the fact that we model this 
variable star at different times. This further supports the requirement of contemporaneous 
observations. Other possible causes for the difference in results could be due the fact that 
the Ha observations are sensitive to different regions within the disk as well as the different 
spatial resolution of the two interferometers. 



There are also other differences between our calculations and iGies et al.l (120071 ) that 
should be highlighted: as noted in Section 2, the calculation of the Ha visibilities and 
the Ha line profiles in our work are obtained from essentially the same calculation that 
determined the radiative equihbrium temperature structure for the disk given the adopted 
density model. The Ha opacity and emissivity at each point in the disk were computed 
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from the local temperature, pressure, and radiative field by solving the set of statistical 
equilibrium equations for a 15-level hydrogen atom. While our fits to the observed Ha line 
profiles are not perfect, they are predicted naturally from the models an d the assurnption 
of pure Keplerian rotation by the disk. To fit the Ha, H7, and Bra lines, iGies et al.l (120071 ) 
found it necessary to convolve their profiles with a unit normalized Lorentzian to broaden the 
line wings; they discuss numerous possible sources for this extra broadening. In the current 
work, we convolved our line profiles only with a unit normalized Gaussian of FWHM of 
0.656 A to bring the resolvin g power of the corn puted profile to 10^ to match the observations. 
Among the reasons cited by iGies et al.l (120071 ) as to possible broad Ha wings, they discuss 
the possibility of Stark broadening for the hydrogen lines, but then note that the densities 
they find are too low for Stark broadening to be important. This is consistent with their 
best fit densities which are about 1-2 orders of magnitudes smaller than those derived here. 

Previous studies have shown that self- consistent disk thermal structures are crucial 



in order to interpre t observations correctly (see for example, ICarciofi fc Bjorkmaru 12006 



Sigut fc Jonesll2007l ). We computed the thermal structure for n Dra with our bedisk code 
us ing the mo d el pa rameters, n = 0.67 ± 0.36, po = 6.2 x 10^^^ gcm~^, R^ = 67R0, listed 



m 



Gies et al.l ( 120071 ). We find substantial variation in disk temperature near the star and 
equatorial plane, and an overall average disk temperature of ~ 13000 K. iGies et al.l (120071 ) 
adopt an isothermal disk t emperature of ~ 9300 K for k Dra b ased on previo u s stu dies 



(ICarciofi fc Bjorkmanll2006l ). but we note that this study as well as lSigut fc Joned (120071 ) do 
not investigate the disk densities as low as the iGies et al.l (120071 ) best-fit model. Finally, we 
note that the value, n = 4.2, we obtain for k, Dra is in agreemen t with typical values previ- 



ously predicted in the literature for Be star disks (see for example IWaterslll986l : iTycner et al. 

mm- 



4.2. p Pisces 



/3 Psc is a B 6Ve star that has lo ng been known to exhibit emission in the hydrogen lines 
(see for example iMerrill et al.lll925l ). This star is similar in spectral type to k, Dra, but it is 
a main sequence dwarf. We have included (3 Psc in this investigation because comparing the 
disk models with the interferometric observations lead to unique density parameters, and it 
is us eful to compare the results fo r this star with k Dra. We adopt the stellar parameters 
from iLevenhagen fc Leister (120041) and note that these parameters are also consistent with 
interpolated values in ICoxl (I2OOOI ). The interferometric observations are from 12 nights 
between 28 Aug. 2005 and 26 Sept. 2005 (see Table H]), and the Ha spectroscopy for P Psc 
(and V Cyg) were obtained on 16 Sept. 2005 and therefore are definitely representative of the 
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Ha line during the period of the interferometric observations. We note that a subsequent 
observation of the Ha hne for j3 Psc on 3 Oct. 2006 had a 20% weaker equivalent width. 

The procedure for selecting the grid of models to be compared in detail with interfer- 
ometric observations was obtained by the same selection process as done for k, Dra. The 
final parameters to be searched covered the range n = 1.5 to n = 5.5 and po = 1.0 x 10~^^ 
g/cm^ to Po = 8.0 X 10^^ g/cm^. In this range, 193 models were constructed for comparison 
with interferometric observations. The corresponding x^ values had a minimum value of 0.83 
corresponding to the model with n = 2.2 and po = 4.0 x 10~^^ g/cm^ and a maximum value 
of 438 corresponding the model with n = 1.7 and po = 1.0 x 10~^^ g/cm'^. The x^ values 
revealed only one other model, n = 2.1 and po = 3.0 x 10^^^, which was statistically the 
same as the best-fit. For this model x^ = 0.84. 

In order to check the predicted Ha profiles, we had to adopt a value for i. The observed 
sharp emission line profile (shown in Figure [6]) suggests a small value for i. However, this 
star has a wsin i of 90° km/s (see Table [2]), so for very small values of i, the star will exceed 
its cri tical velocity . The estimated critical velocity of a star of this spectral type is ~ 416 



km/s (jPorteiill996l ). Therefore, we argue that the inclination of this system must be at least 



15°, and we have adopted a value of i = 20°. 

We calculated the Ha equivalent widths for a subset (~ 20 models) chosen by the 
best-fit x^ values, and we considered x^ values up to ~ 1.52. Statistically, the two models 
corresponding to x^ values of 0.83, and 0.84 are equally as good while the other models in the 
subset with larger values of x^ are poorer fits to the interferometric data. Nevertheless, to be 
complete we computed profiles for every model in the subset. The equivalent widths for the 
two models corresponding to the minimum x^ values had reasonable line shapes compared 
to the observed profile but the model corresponding to n = 2.2 had an equivalent width of 
-25.0 A while the model corresponding to n = 2.1 had an equivalent width of -21.6 A which 
is closer to the value for the observed spectra of -15.3 A. 

Figure [6] shows the Ha profile for the best-fit model {n = 2.1, po = 3.0 x 10~^^ g/cm^) 
in the left-most panel in terms of interferometric observations and spectroscopy. Models 
corresponding to other values of n did not produce Ha profiles that agreed with observations. 
Figure [6] shows the behavior of the Ha profile with increasing values of n. The reader can 
clearly see that the wings of the line in the middle and right-most panels are too wide, 
resulting in predicted emission which is too large. These profiles correspond to models with 
n = 3.6, Po = 1.5 X 10"^° g/cm^, equivalent width of -30.4 A, and n = 4.0, po = 4.0 x 10"^° 
g/cm^, equivalent width of -23.4 A, and have interferometric reduced x^ values of 1.13, 
and 1.21, respectively. Models with values of n < 2.1 had large x^ values based on the 
interferometric results, and were eliminated. Figure [7] shows the interferometric observations 
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Fig. 6. — The observed Ha line for (3 Psc compared with model predictions for a series of 
models showing the changes in the profile as n is increased. The model profiles are the solid 
lines and the circles represent the observation. 



-19- 



for our best fit model with n = 2.1, and po = 3 x 10 ^^ g/cm^ for j3 Psc. 



4.3. V Cygnus 



V Cyg has long been recognized as a Be star with strong Ha emission ( lFleming)ll89ll : 
Campbell! Il895l ). Although several authors have noted th at this emiss ion stays relatively 
constant over periods as long as decades (see, for example, |Peterslll9791) others have noted 
perio dic outbursts in Ha emission, typical of many Be stars (JBallereau et al.l Il995l : iNeiner 
20051 ). The interferometric observations for v Cygnus are from 12 nights between 28 Aug. 
2005 and 26 Sept. 2005, and the Ha line profile was acquired on 16 Sept. 2005 (see Tabled]). 
We note that we also have observations from 17 Oct. 2005 that confirm the spectroscopic 
stabihty of the star. 



Neinen ( l2005l ) suggested that the outburst s, which occu r on periods of ~ years, are the 



result of multi-periodic, non-radial pulsations. iNeinerl (120051 ) derived the stellar parameters 
for this star by detailed modeling constrained by photometric and spectroscopic data. Three 
of their models include veiling effects, and gravitational darkenin g for a range of stellar 
rotational velocities, 0.80, 0.90, and 0.95 of the critical velocity. (See lNeinerl (120051 ) for more 
detail.) We have chosen to adopt the stellar parameters from their model D which has an 
intermediate value of the rotational velocity (that is, 0.90 of critical velocity). Also. iNeiner 
(120051 ) finds a value for the inclination angle of i ~ 27 °. The stellar radius listed in Table [2] 



corresponds to the equatorial radius given by JNeineri (120051 ). This is a reasonable choice 
since the star is nearly pole-on. The stellar parameters adopted for v Cyg are reasonable 
conipared to other values presented in the literature of this spectral type (see for example 
Co2\200d \. 



The procedure for selecting the grid of models to be compared with interferometric 
observations was obtained by the same selection process as done for the two other stars in 
this investigation. The final selection of parameters to be searched ranged from n = 1.8 
to n = 4.5, and po = 1.0 x 10~^^ g/cm^ to po = 4.0 x 10~^g/cm^, giving a total of 193 
models constructed for comparison with interferometric observations. The corresponding 
X^ values had a minimum value of 1.14 corresponding to the model with n = 4.0 and 
Po = 8.0 X 10^^° g/cm'^ and a maximum value of 471 corresponding the model with n = 1.8 
and Po = 8.0 x 10^^^ g/cm^. The x^ values revealed several models that statistically fit the 
observations equally well. Again, we turn to spectroscopy to remove the degeneracy. 



In order to predict Ha profiles we adopt i = 30° which is consistent with JNeinerl ( 120051 ) . 
We calculate the Ha equivalent widths for a subset (~ 22 models) chosen based on the 
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Fig. 7. — The interferometric data and model with n = 2.1, and po = 3.0 x 10~^^ g/cm^ for 
(3 Psc. The dashed hne represents the central star and the solid line is obtained by taking a 
Fourier transform of a synthetic image of the preferred model. 
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interferometric reduced x^ < 2. Statistically, the models with larger values of x^ are poorer 
fits to the interferometric data. The equivalent widths for the subset ranged from a minimum 
of -19.0 A to a maximum of -43.1 A. The observed Ha line has an equivalent width of -24.8 
A, and there were three models which had reasonable line shapes and equivalent widths 
that agreed within ±2 A and had statistically similar interferometric ~x^ values. Figure [8] 
shows the Ha profiles corresponding to these models. Table [3] lists the values of n, po, Ha 
equivalent width, and the reduced x^ values based on comparisons from interferometry for 
these 3 models. The models with n = 4.1 and n = 4.2 in Figure [H] have predicted profiles 
that are too narrow and also do not fit the peak of the observed profile. Since the x^ values 
are statistically similar, we adopt the model corresponding to n = 4.0 and po = 8.0 x 10^^'^ 
g/cm^ in Figure [8] as the best-fit. As a further illustration. Figure [9] shows Ha profiles 
corresponding to models with smaller and larger values of n than those given in Table [3l 
Specifically, these profiles correspond to models with n = 3.6, po = 3.0 x 10~^° g/cm^ and 
n = 4.4, Po = 2.0 X 10~^ g/cm^ with reduced x^ values of 1.37, and 1.29, respectively. The 
shapes of these profiles do not match the observed line. The predicted profile corresponding 
to n = 3.6 has too much emission (with an equivalent width of -33.7 A) and is too wide in 
the line wings. The predicted profile corresponding to n = 4.4, has too little emission (with 
an equivalent width of -20.3 A) and much of the profile is narrower than the observation. 
Figure [10] shows a comparison of the interferometric observations and the best-fit model for 
V Cyg. 



5. Conclusions 

In this work, disk density models were obtained for the classical Be stars k Dra, (3 Psc, 
and V Cyg by matching the observed interferometric Ha visibilities with Fourier transforms 

Table 3: Best-Fit Models for v Cygnus 



n 


Po 


Ha Equivalent 


^2 a 




g/cm^ 


Width [A 




4.0 


8.0 X 10-1° 


-26.0 


1.14 


4.1 


1.0 X 10^9 


-24.4 


1.16 


4.2 


1.5 X 10-9 


-23.7 


1.30 



"Based on comparison with interferometry. 
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Fig. 8. — The observed Ha line for v Cyg compared with the best predicted hnes from 
the models which are also in agreement with the interferometric observations. The model 
profiles are the solid lines and the circles represent the observed line. See Table |3] for the 
model parameters, Ha equivalent widths, and the reduced x^ values from interferometry 
corresponding to these models. 
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Fig. 9. — The observed Ha lines for v Cyg compared with the best predicted hnes for n = 3.6 
and n = 4.4. The model profiles are the solid lines and the circles represent the observed 
line. The model on the left panel has a lower value of n and higher density, resulting in an 
Ha profile which is too large. The model on the right panel has a higher value of n and 
lower density, resulting in an Ha profile which is too small. 
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Fig. 10. — The interferometric data and model with n = 4.0, and po = 8.0 x 10^^° g/cni'^ for 
V Cyg. The dashed hne represents the central star and the solid line is obtained by taking 
a Fourier transform of a synthetic image of the preferred model. 



-25- 



of synthetic images produced by theoretical models. It was demonstrated that the additional 
constraint obtained by the (co-temporal) observed Ha line profile was critical to selecting 
a density model. This was particularly evident in the case of k Dra. We have demon- 
strated that in the case of these 3 stars, only by additional observational constraints can the 
degeneracy in the interferometric best-fitting models be removed. 

The best fit models to both the interferometric observations and the observed Ha line 
profiles for the three stars are summarized in Table HI The range of n is consistent wi th 



previous determinations based on other diagnostics (such as IR excess see - IWaterd ( 1l986l )). 
We note that the range of "base" densities, po, varies by over two orders of magnitude 
between the three stars. 
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Laboratory and the US Naval Observatory, in cooperation with Lowell Observatory, and 
is funded by the Office of Naval Research and the Oceanographer of the Navy. We thank 
the Lowell Observatory for the telescope time used to obtain the Ha line spectra presented 
in this paper. This research was supported in part by NSERC, the Natural Sciences and 
Engineering Research Council of Canada. We thank L. Thomson for contributions to this 
work as an NSERC USRA in 2006. Facilities: LO:42in(), NPOI () 
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